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ABSTRACT 

We present the results of a uniform analysis of the soft X-ray spectra of fifteen type I 
AGN observed with the high resolution X-ray gratings on board Chandra. We found 
that ten of the fifteen AGN exhibit signatures of an intrinsic ionized absorber. The 
absorbers are photoionized and outflowing, with velocities in the range 10^ — 10'^ km 
s~^. The column density of the warm absorbing gas is ~ 10^"^^'^ cm^^. Nine of the ten 
AGN exhibiting warm absorption are best-fit by multiple ionization components and 
three of the ten AGN require multiple kinematic components. The warm absorbing gas 
in our AGN sample has a wide range of ionization parameter, spanning roughly four 
orders of magnitude ~ 10'^"^ ergs cm s~^) in total, and often spanning three orders 
of magnitude in the same gas. Warm absorber components with ionization parameter 
^ < 10 generate an unresolved transition array due to Fe in seven of the ten AGN 
exhibiting warm absorption. These low ionization state absorbers may also carry away 
the largest mass outflows from the AGN. The mass outflow rate depends critically on 
the volume filling factor of the gas, which cannot yet be directly measured. However, 
upper limits on the mass outfiow rates for filling factors of unity can be much greater 
than the expected accretion rate onto the central supermassive black hole and filling 
factors as small as 1% can give outflow rates comparable to the accretion rate. There 
appears to be a gap in the outflow velocities in our sample between ~ 300 — 500 km 
s~^, the origin of which is not clear. The outflow components with velocities below this 
gap tend to be associated with lower column densities than those with with velocities 
above the gap. 

Key words: galaxies: active - galaxies: individual - galaxies: Seyfert - techniques: 
spectroscopic - X-rays: line - emission: accretion - disks :galaxies 



1 INTRODUCTION 

X-ray emission from active galactic nuclei (AGN) is believed 
to be powered by an accretion flow onto a supermassive black 
hole (SBH). In the unified AGN paradigm, type I AGN have 
accretion disks at a small angle of inclination to the ob- 
servers' line-of-sight. Hard X-ray (>2 keV) spectra from 
type I AGN typically reveal a continuum that is well de- 
scribed by a simple cut-ofi' power-law model, with a fluo- 
rescent Fe K line complex (and sometimes a reflection con - 
tinuum) superimposed (see eg. [Reynolds fc NowakI (|2003l ) 
and references therein). The soft X-ray (<2 keV) spectra of 
type I AGN are typically complex and can yield a great deal 
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of information about the distribution and state of matter in 
the AGN central engine. 

Spectral complexity beyond a simple power law model 
in the soft X-ray spectrum of an AGN was first observed 
in MR 2251-178 iHalpernl (Il984l '). The spectral complex- 
ity was due to absorption by partially ionized, optically 
thin, circumnuclear material. This material, termed the 
'warm absorber' was proposed as a common constituent of 
many AGN. ROSAT observations of nearby type I AGN 
show ed that warm absorption was not uncomm on, see 
e.g. (|Nandra fc Founds! Il99^ : [Turner etHI 1 19931 ). Subse- 
quent observations with ASCA detected O vii and O viii 
absorption e dges due t o wa rm absorption in ~ 50% of 
type I AGN (|Revnoldsl (|l997D . George et al. 1998). Obser- 
vations in the UV band reveal warm absorption with mul- 
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Table 1. The Chandra HETGS sample of type I AGN . AGN are listed in order of increasing RA. Columns 2,3 and 4 
give the co-ordinates and the redshift of the source (from NED). Redshift was deduced from observations of the 21cm H I 
line where possible, since optical estimates of z may be confused by AGN outflo w. Column 7 is th e mean count rate for 
the combined ±1 ord er spectra of the Medium Energy Grating (MEG). " from lElvis et all lll989l ). except for Mkn 509 
iMurphv et al] lll996l '). The Galactic column dens ity toward s F9, NGC 7314, NGC 3516, NGC 3783, NGC 5548, Mkn 
766, NGC 3227 and Akn 564 was obtained from lRevnoid3 lll997l'l . NGC 3516: Dates of three 'snapshots' are given; 
NGC 3783: data were combined from five observations over a period of 124 days; MCG— 6-30-15: this observation was 
made over a period of ~ 138 days in three parts and only the combined data are analyzed here. Total good integration 
time of spectrum. 3C 120 is also classified as a broad— line radio galaxy (NED). To gauge changes in AGN spectral 
continuua, particularly in the vicinity of the Oxygen absorption edges, the XMM-Newton EPIC PN spectrum of this 
AGN was studied. ^ we summed the first two observations of NGC 3516 into a single 'low state' spectrum and compared 
this with the 'high state' spectrum in the third observation. 
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tiple velocity components and ionization states, although 
the relationship between the UV absorbers and the X- 
ray absorbers is s till not clear l|Crenshavy fc Kraemer|[l999l : 
lAravet al.ll2002h . The unprecedented spectral resolution 
of the gratings on board Chandra and XMM-Newton has 
permitted the detection of discrete soft X-ray absorp- 
tion features and emission lines for the first time. The 
resulting picture of the X-ray warm absorber in many 
type I AGN is of an outflow exhibiting multiple narrow 
absorption linos correspond ing to di fferent ionization states 
e.g. ( ColhuKc ct al. 2001; Lee et al.1 [20 01: Sako ct al. 2001; 
Kaastra et all 120021: INetzer et al.ll2003l : lYagoob et al.l 12003 : 
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The present generation of X-ray 
detectors also have the spectral resolution to detect unre- 
solved transition arrays in moderately ionized Iron (Fe""^^) 
in several AGN spe ctra (Sako ct al. 2001: Behar et al.ll200il: 
iBlustin et allbooi iPounds et al.ll2003h . 

Our picture of the X-ray warm absorber is now more 
complicated, and many fundamental questions remain. Does 
the warm absorber 'kn ow' about the SBH mass, e.g. 
[Morales fc FabianI l|2002l ), the accretion rate, or the AGN 
luminosity? Does radiation pressure dominate the warm ab- 
sorber outflow? Is there a link between the warm absorber 
and the Fe K band emission, e.g. Matt et al. (1994)? Only 
via the analysis of the soft X-ray spectra of a sample of 
type I AGN can we answer such questions. Here we present 
the results of a uniform analysis of the soft X-ray data from 
a sample of fifteen type I AGN observed with the high en- 
ergy transmission grating spectrometer (HETGS ) on board 



Chandra (|Markert et al.ll 19951 ). Our aim is to study the soft 
X-ray spectra of these AGN so that we may begin to an- 
swer some of the outstanding questions about warm ab- 
sorption in AGN. Of course our study is limited by both 
our choice of sample and X-ray instrument. The Chandra 
HETGS bandpass is less sensitive to low ionization state 
absorbers than the Chandra LETGS bandpass for exam- 
ple. Both the LETGS and the RGS (aboard XMM-Newton) 
have a higher effective area than the MEG at low energies 
but the HETCS spectral resolution is superior to both the 
LETCS and RGS. Utilizing the best spectral resolution cur- 
rently available is the prinicipal driving factor for using the 
HETGS in the present study. There are many more HETGS 
observations than LETCS observations and the RGS band- 
pass does not extend to the Fe K region of AGN spectra. We 
note that in a uniform analysis, individual source pecular- 
ities (e.g. in the continuum modelling) may be missed. An 
additional complication is that the absorbers in some AGN 
are known to vary, so drawing general conclusions based on 
a snapshot of a variable absorber may not be warranted in 
some cases. Nevertheless, a uniform analysis is useful since 
variations between analysis software and/or methodologies 
can account for considerable differences in interpretation, 
often over the same data. 



2 THE SAMPLE AND DATA ANALYSIS 

Our study is based on the sample of fifte en type I AGN 
selected bv lYagoob fc PadmanabhanI (|2004l ). The AGN are 
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listed in Table [T] were originally assembled for a study of the 
Fe K band emission and had z < 0.05, a total, first-order, 
high energy grating count -rate of > 0.05 ct/s and the ob- 
servations were in the Chandra public data archive as of 
July 1, 2003. This constitutes a rather heterogeneous sample 
that is not based on a scientifically motivated selection cri- 
terion. Therefore, certain results and conclusions pertaining 
to sample properties (such as the fraction of AGN exhibit- 
ing signatures of photoionized outfiows) must be interpreted 
with the appropriate caution. However, members of the sam- 
ple do satisfy one very important criterion relevant for this 
study, namely that these sources are some of the bright- 
est members of their class and therefore lend themselves to 
performing detailed X-ray spectroscopy with the Chandra 
HETGS. This is not a coincidence because generally speak- 
ing, the highest signal-to-noise members of a class tend to 
get accepted first by selection panels for observations in the 
early years of a new mission. 

We note that several of the AGN listed in Table [T] have 
been observed with Chandra again (for which that data 
became public after July 1, 2003). Furthermore, an addi- 
tional thirteen AGN have been observed with the Chan- 
dra gratings that would fulfil l the selection criteria of 
lYaaoob fc Padmanabhaiil (|2004l '). In the future we intend to 
extend our sample study to include these AGN & more re- 
cent observations of the AGN in our sample. Also listed in 
Tableware the AGN redshifts (from NElfl using 21cm H i 
radiation measurements where possible), the R A and DEC 
(also from NED), the Gal actic column density (|Elvis et al.l 
1 19891 : iMurphv et al.|[l99^ and the total exposure times of 
the spectra. 

The Chandra data were reprocessed using ciao 2.1.3 
and CALDB version 2.7, according to recipes described in ciao 
2.1.3 thread^. The instrument in the focal plane of Chan- 
dra during the observations was the HETGS, which consists 
of two grating assemblies, a high-energy grating (HEG) and 
a medium-energy grating (MEG). Only the summed, nega- 
tive and positive, first-order Chandra grating spectra were 
used in our analysis. The HEG bandpass is ~0.8-10 keV and 
the MEG bandpass is ~0.5-10 keV but the effective area of 
both instruments falls off rapidly at either end of the band- 
pass. Since the MEG soft X-ray response is much better than 
the HEG we used the MEG as the primary instrument. 

We made effective area files (ARFs or ancillary re- 
sponse files), p hoton spectra and c ounts spectra following 
the method of lYagoob et al.l (|2003l ') . We did not subtract 
detector or X-ray background since it is such a small frac- 
tion of the observed counts. For spectra with zero or few 
counts per bin anywhere in the bandpass, attempting to 
subtract background whilst retaining the best spectral res- 
olution possible can result in worse systematic errors com- 
pared to the case when no background subtraction is at- 
tempted. This is especially true when the background itself 
is weak, having zero counts for most spectral bins. However, 
background could be a source of contamination at the low- 
est energies of the MEG spectra (where the effective area 
is smallest), for weak and/or heavily absorbed sources. We 
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examined spectra taken from two strips, either side of the 
on-source data, to check the level of background for each 
data set. We found that contamination could be a problem 
in NGC 3516, NGC 3227 and NGC 7314, where the back- 
ground level becomes comparable to the source intensity be- 
low ~ 0.6 keV. We will bear this in mind when interpreting 
the data. We note that five of the fifteen AGN in our sam- 
ple were observed with the Chandra low energy transmission 
grating spectrometer (LETCS). The high resolution camera 
(HRC) was used in four of these five LETCS observations 
(as opposed to ACIS). However, order separation is not pos- 
sible with the HRC and the properties of the LETCS are 
considerably different to those of the HETGS (e.g. LETCS 
has a lower spectral resolution of only 0.05A FWHM ) and 
so inclusion of the LETCS data would have complicated ef- 
forts to perform a uniform analysis and comparison amongst 
the sources. Note that the observations of three of the AGN 
(NGC 3516, NGC 3783 and MCG -6-30-15) were not made 
in a single observing period. NGC 3516 was observed in three 
parts, as detailed in Table [T] However, two of the observa- 
tions found NGC 3516 in a low fiux state. So, in our analysis 
of the warm absorber in this AGN, we summed two of the 
observations into a single 'low state' spectrum and compared 
this with the 'high state' spectrum in the third observation. 
NGC 3783 was observed in five snapshots over a period of 
~ 124 days. However, we analyze only the summed data in 
the present paper, since t he spectral v a riabil ity, which has 
been studied in detail bv iNetzer et al.l l|2003l ) is not great. 
The observation of MCG —6-30-15 has a single sequence 
number and observation ID, but was made in three parts, 
spanning a period of ~ 138 days. The S/N in the individual 
parts was not sufficient to warrant analyzing the three parts 
separately so here we analyze only the summed data. 

For comparison with some of the results from the Chan- 
dra observations, we studied XMM-Newton spectra from 
some of the AGN in our sample. We used observations made 
with the XMM-Newton EPIC FN instrument. The data used 
were obtained from the reduced data products from the Xas- 
sist database 0. AGN in our HETGS sample with corre- 
sponding XMM-Newton observations are listed in Table [1] 

We treated the statistical errors on both the pho- 
ton and counts spectra with particular care since the low- 
est and highest energies of interest can be in the Pois- 
son regime, with spectral bins often containing a few, or 
even zero counts. When plotting data, we assign statis- 
tical upper and lower errors of 1.0 -\- ^/{N-\-0. 75) and 
A^(1. 0- [1.0- (l/(9Af)) - respectively (jCehrelsl 

Il986h on the number of photons, A'', in a given spectral 
bin. When fitting the Chandra data, we used the C-statistic 
(|Cashlll976 ') for finding the best-fitting model parameters, 
and quote 90% confidence, one-parameter statistical errors 
unless otherwise stated. The C-statistic minimization algo- 
rithm is inherently Poissonian and so makes no use of the 
errors on the counts in the spectral bins described above. 
All model parameters will be referred to the source frame, 
unless otherwise noted. Note that since all models were fit- 
ted by first folding through the instrument response before 
comparing to the data, the derived model parameters do not 
need to be corrected for instrumental response. 

* |http : //xassist .pha. jhu. edu| 
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We used XSPEC v. 11. 3.1 for spectral fitting to the 
HETGS spectra. All spectral fitting was done in the 0.5- 
5 keV energy band, excluding the 2.0-2.5 keV region, which 
suffers from systematics as large as ~ 20% in the effec- 
tive area due to limitations in the calibration of the X-ray 
telescope 0. We performed spectral fits using data binned 
at ~ 0.02A(which is close to the MEG FWHM spectral 
resolution of 0.023A), unless otherwise stated. There has 
been a continuous degradation of the quantum efficiency 
of Chandra ACIS with time, due to m olecular contamina- 
tion B In anal yses of individual sources (|Yaaoob et al.ll2003l : 
iMcKernan eF al. 2003a bj we found that a pure ACIS correc- 
tion (the worst-case effect) affects only the inferred intrinsic 
continuum (at less than ~ 0.7 keV) and does not affect the 
important physical parameters of warm absorber models. A 
detailed comparison of spectral results with and wit hout the 
ACIS degradation correction was given in (McKcrn an et all 
1200381 ) for the case of NGC 4593. We also note that the data 
in our sample are from early in the Chandra mission when 
the ACIS degradation was not so severe and the effects of 
the degradation are mitigated by the limited signal-to-noise 
ratio of the data at the lowest energies. 



3 THE OVERALL SOFT X-RAY SPECTRA 

The soft X-ray continuua of Type I AGN can be quite com- 
plex. An excess over a simple power-law continuum model, 
as well as possible absorption features, are common. In 
our sample, nine of the fifteen AGN in Table [T] show fea- 
tures characteristic of a soft excess and/or absorption edges. 
NGC 3227, NGC 3516, IC 4329A and NGC 7314 are heavily 
absorbed in the soft X-ray band. F9, 3C 120 and Mkn 279 
have relatively simple continuua, since their spectra show 
little or no evidence for a soft excess. Since the soft excess 
typically appears only in the 0.5-0.7 keV band of our data, 
we do not have enough information to constrain its origin 
and so sophisticated modelling of the continuum (such as 
with a power-law plus a blackbody model component) is 
not warranted. We found that a broken power-law (although 
likely to be non-physical) is adequate to describe the intrin- 
sic continuum including any soft excess in the HETCS en- 
ergy band. When partially covering cold absorption is added 
to the model a useful empirical description of the overall 
spectra is obtained, leaving residuals mostly due to complex, 
ionized absorption, and in some cases some emission lines. 
Table 2 shows the best-fitting parameters obtained with this 
model for each source, as well as observed 0.5-2 keV fiuxes. 



4 PHOTOIONIZATION MODELING OF THE 
AGN WARM ABSORBERS 

We used the publicly available photoionization code XS- 
TAR 2.1.kn3[j to generate several grids of models of emis- 
sion and absorption from photoionized gas in order to di- 
rectly compare with the data. Version 2.1.kn3 includes un- 
resolved transition arrays (UTAs) of inner-shell transitions 



of Fe for the first time in XSTAR. It is important to include 
UTAs in our photoionized model grids since they provide the 
key to studying low ionization state absorbers in AGN and 
UTAs of moderately ionized Fe""^^^ have bee n observed 
in th e spectra of several t ype I AGN ( see e . g. llSako et al 
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[2003)). We use d the default solar ab undances in XSTAR (e.g. 
see Table 2 in lYa aoob et all |2003)). Note that the XSTAR 
line database has been constructed from lines withpublished 
wavelengths, such as the Chianti line com pilation [f|. This al- 
most completely excludes the lines listed in lBehar fc Neztzerl 
(2002) for example, which are based on HULLAC calcula- 
tions and which are only reliable at about the 1% level. 
These lines, due to L-shell transitions in Ne, Mg, Al, Si, S, 
Ar, Ca a nd Fe, are pre s ent in many Chandra grating spectra 
(see e.g. iNetzer et al.l (|2003l )). but so far they have mostly 
been omitted from the XSTAR data base rather than com- 
promise the accuracy of the database iKaUmanI (j2003l ). 

There may be additional absorption in most of these 
AGN, beyond the simple photoionized warm absorber mod- 
els that we consider in this uniform analysis. Candidates 
for the extra absorption include neutral dust ( e.g. possi- 
bly in the form of Fe02 in MCG-6-30-15, see iLee et all 
(I2OOII )') or deeper O Vll or O vill edges due to Oxygen over- 
abundance. Alternatively, a more complex model of the con- 
tinuum may be more appropriate, incorporating relativisti- 
cally broadened emission lines in the soft X-ray band (see 
e.g. (iBranduardi-Ravmont et all I2001I : iMason et all l2003l : 
ISako et al.ll2003l )V However, investigation of these issues is 
beyond the scope of the present uniform analysis and we 
shall return to these important points in future work. 

The best spectral resolution of the MEG cannot directly 
constrain turbulent velocities with a "b-value" of less than 
170 km s~^. Furthermore, the XSTAR model spectra are cal- 
culated by XSTAR on a grid that does not always preserve 
accuracy in the line widths and line equivalent widths (Kail- 
man, private communication). Therefore, direct fits to the 
data using the XSTAR spectra are approximate in the re- 
gions containing absorption lines. However, the ionic column 
densities output by the XSTAR are more accurate than the 
XSTAR spectra and can be used for more detailed modeling 
of individual absorption lines. This level of detail is beyond 
the scope of the present uniform analysis. We used a velocity 
turbulence (b-value) of 170 km s~^, which corresponds ap- 
proximately to the limiting MEG spectral resolution (~ 300 
km s^^ FWHM at 0.5 keV), given that the data cannot 
directly constrain smaller line widths. It is possible to con- 
strain b using a curve-of-growth analysis if one has several 
absorption line measurements, and this has been done for 
some of the sources in our sample in more detailed stud- 
ies. For example, in Mkn 509 b is consistent with ~ 100 km 
s"^ ijYaaoob et al.l (|2003l )). with similar results obtained for 
other sources. Since 6 is generally comparable to or less than 
the spectral resolution, our adopted value of b is justified for 
the purpose of spectral fitting. 
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Figure 1. Spectral energy distributions (SEDs) used for plio- 
toionization modelling to the Chandra data of ten AGN in our 
sample (see i|4.1| l. In order of decreasing luminosity at lOkeV 
(right-hand-side of graph), the SEDs correspond to Mkn 509, IC 
4329A, NGC 5548, NGC 3783, Mkn 766, Akn 564, NGC 4593, 
NGC 3516, MCG-6-30-15 and NGC 4051 respectively 



4.1 XSTAR modelling procedure 

First, we constructed a spectral energy distribution (SED) 
for each A GN in our sample, a ccordiuK to the method de- 
scribed in lYaaoob et all l|2003h . In principle, derivation of 
the SED should be an iterative process since the intrinsic 
X-ray spectrum is obtained from fitting the data using an 
SED that must already contain information about the X-ray 
spectrum. In practi ce, the method employed (described in 
lYaaoob et all ()2003l )) is a two-step process in which the first 
step is an initial estimate of the intrinsic X-ray spectrum 
is obtained by fitting the X-ray data with a power-law or 
broken power-law only and no photoionized absorber. The 
resulting estimated SED is used to generate grids of pho- 
toionization models that are used to obtain better fits to 
the data and to derive a new intrinsic X-ray spectrum. This 
is then used to generate another SED which is in turn used 
to generate a new set of photoionzation model grids and 
the latter are the ones used for the final model fits. Aver- 
age radio, infrared (IR), optical and u ltraviolet (UV) fluxes 
were obtained from lward et al.1 (Il987l ') and NED. Note that 
flux measureme nts, at differ e nt wa velengths, obtained from 
historical data ijWard et al.l l|l987i ). NED) are not in gen- 
eral contemporaneous. For each source, the 0.5 keV intrin- 
sic model flux from the intrinsic X-ray continuum model 
was then simply joined onto the last point of the UV part of 
the SED by a straight line in log-log space. The hard X-ray 
power law was extended out to 500 keV. Generally, the high- 
energy cut-off in the observed X-ray spectra of Seyfert galax- 
ies lies in the range 100-500 keV; only blazars and BL Lacs 
are observed to have significant flux beyond 500 keV. The 
ionization balance is is not sensitive to the exact position of 



the cut-off between 100-500 keV . At the low-energy end, it 
has been shown (for example, bv lFerland et al.l (|2002l )) that 
if the observed IR continuum carries a signiflcant thermal 
component from reprocessing of the instrinsic continuum, 
using the observed SED as an input to photoionization mod- 
els can potentially affect the resulting ionization balance of 
the plasma. However, our model fits are driven by X-ray 
features and we showed in a detailed study of NGC 4593 
ijMcKernan et all (l2003al ')') that removing the prominent IR- 
optical continuum bump from the SED yielded warm ab- 
sorber parameters that were within the 90% confidence in- 
tervals obtained when the bump was not removed. Figure [T] 
shows the SEDs that we used for the ten AGN in our sample 
that exhibited warm absorption (see 814.21 below). 

The photoionization model grids used here are two- 
dimensional, corresponding to a range in values of total neu- 
tral Hydrogen column density, Nu, and the ionization pa- 
rameter, ^ = Lion/irieR^)- Here Lion is the ionizing luminos- 
ity in the range 1-1000 Rydbergs, Ue is the electron density 
and R is the distance of the illuminated gas from the ion- 
izing source. The ionizing luminosities were calculated from 
the SEDs by normalizing the absorption-corrected 0.5-2 keV 
fluxes. The grids were computed for equi-spaced intervals 
in the logarithms of A^h and log^, in the ranges 10^^ to 
5 X lO^^cm"^ and —1.0 to -f4.0 erg cm s~^ respectively. We 
computed grids with Tie in the range lO^"^'^ cm~^. For XS- 
TAR models of the absorber, we conflrmed that results from 
fitting the photoionization models to the X-ray data were 
indistinguishable for densities in the range Ue — 10^ cm~'^ 
to 10^^ cm~^ for all of the AGN. Hereafter we will use 
Tie = 10* cm~^ unless otherwise stated. 

Our aim was to fit the MEG spectra with a simple 
power-law or a broken power law continuum (whichever was 
the better fit), modified by absorption from photoionized gas 
(derived from the AGN SED), neutral gas intrinsic to the 
AGN and Galactic absorption. All model fitting to the data 
was carried out using XSPEC vll.3.1. In order to fit XS- 
TAR photoionization model grids to the data, we required 
an offset velocity for the warm absorber for each of the AGN. 
Therefore, first we tested for the following He-like and H-like 
absorption lines in each spectrum within ~ ±5000 km s~^ of 
the line rest-energy in the AGN frame: N vil Lya (A24.781 
A), O VII (r) (A21.602 A), O viii Lya (A18.969 A), Ne ix (r) 
(A13.447 A), Ne x Lya (A12.134 A), Mg xi (r) (A9.169 A), 
Mg XII Lya (A8.421 A), Si xill (r) (A6.648 A), Si xiv Lya 
(A6.182 A). Note that for individual absorption features, the 
statistical errors on the line energies are typically very small 
and so the systematic uncertainty likely dominates these er- 
rors. The relative wavelength accuracy of the MEG is 0.0055 
A and the absolute wavelength accuracy of the MEG is 0.011 

m 

Next, we established an offset velocity for a particular 
AGN by finding the weighted-mean offset velocity (and the 
90% confidence limits) of the centroids of all the absorption 
lines detected at > 90% significance in the spectrum. Some 
AGN exhibited groups of two or more absorption features 
separated from another such group by more than the 90% 
errors on their respective centroid velocities. In these cases, 
we established more than one weighted-mean offset velocity. 



3 http://asc.harvard.edu/proposer/POG/html/HETG.html 
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Finally, we fitted XSTAR models only to those AGN where 
two or more prominent absorption signatures had similar off- 
set velocities (within 90% confidence limits). This approach 
has the merit of being uniform, but it runs the risk of be- 
ing too conservative and we may be biasing our study in 
favour of more prominent warm absorption signatures. On 
the other hand, our method does have the merit of ignoring 
statistically spurious features. We tested for the presence of 
multiple ionization components in the absorbers by investi- 
gating the statistical significance of additional model compo- 
nents within the 90% confidence limits of the weighted-mean 
offset velocity. 



4.2 XSTAR model results 

Table 13] lists the best-fit warm absorber models that we es- 
tablished for our AGN spectra. Figures [2131 and l4l show the 
best-fit models overlaid on the HETGS spectra. We note 
that in Table [3] the derived column densities of each warm 
absorber component do not depend on the "b-value" as- 
sumed in the XSTAR models because the absorption lines 
are unresolved. Also, these column densities are not soley 
driven by the absorption lines. In particular, in the case of 
NGC 3783, altough our simplistic model gives a poor fit in 
some regions of the spectrum, missing some absorption lines 
altogether, the column densities of the four warm absorber 
components are not sensitive to this since the four numbers 
are determined by a fit to tens of absorption lines as well as 
to the overall spectral shape. 

Only ten of the fifteen AGN in our sample satisfied our 
criteria for selection and spectral fitting with XSTAR grids. 
The AGN F9, 3C 120, Mkn 279, NGC 3227 and NGC 7314 
are not included in Table [3] In the case of F9, we found 
weak O VIII Lya and O Vll (r) absorption signatures 90% 
confidence). 3C 120 exhibited strong O vill Lya absorption 
but no other absorption signatures at > 90% confidence. In 
Mkn 279, Ne x Lya was the only absorption feature present 
at > 90% confidence. In NGC 3227 and NGC 7314, the 
soft X-ray band was very strongly absorbed, so we could 
not distinguish between a cold absorber and a mildly warm 
absorber. In NGC 3227 we only detected a single statisti- 
cally significant absorption feature (Si xiv Lya at ~ —1700 
km s~^ relative to systemtic velocity). Note that previous, 
higher S/N data of the heavily absorbed AGN in our sample 
(NGC 3227 & NGC 7314) have in some c ases detec t ed ev i- 
dence for photoio ni zed absorption (see e.gl Revnoldj l|l997h : 
iPtak et al.1 (|l994l ): ICeorge et al.1 (|l998l )). Spectral plots of 
the data for the five AGN that are not shown in Figure [3] 
can be found in the HETG public database, HotGA^ 

Table U shows results from the literature for model fits 
to X-ray spectra from the AGN in our sample. In general, 
given the limitations of our methodology, the agreement is 
actually quite good. Our uniform analysis ignores individual 
continuum peculiarities in individual AGN and our selection 
criteria for absorbers are more conservative than many in the 
literature. In spite of this, our method seems to pick up the 
same significant absorber components in the literature. Per- 
haps the best example of this can be seen by comparing our 
results from NGC 3783 in Table [3] with those from a detailed 



' |http : //hotgas .pha. jhu. edu| 



analysis bv lNetzer et al.l (l2003l) of the best S /N X -ray spec- 
trum yet. Table l4l shows that iNetzer et al.l l|2003l ') find two 
velocity components of the warm absorber with three ioniza- 
tion states. Our simple, uniform method picks out the same 
velocity components at ~ —500, —1100 km s~^ and reveals 
the three ionization states, although we do not establish 
high velocity components for the middle and low ionization 
states. The total column density of the absorbers a grees very 
well w ith the value of ~ 3 x lO^'^cm"'^ found bv lNetzer et al.l 
(|2003h . Our simple, uniform analysis m eans that our contin- 
uum fit is not as good as that found bylNetiereFaD (|2003h . 
but it is encouraging that we reproduce most of their model 
parameters, given that their analysis of this spectrum was 
far more detailed than ours. 

Of course, there are AGN for which individual source 
peculiarities mean that our results agree less well with those 
in the literature. In the case of MCG -6-30-15 . iLee et al.1 
(|200ll ') and lSako et al.l l|2003l ) come to very different conclu- 
sions concerning both the origin of the continuum in this 
AGN and the parameters associated with the warm ab- 
sorber. 

Unresolved transition arrays (UTAs) of Fe are a key di- 
agnostic of material in a low-ionization state, which could 
account for the bulk of the mass outflow from AGN (depend- 
ing on the absorber density and geometry). The present gen- 
eration of X-ray telescopes have the spectral resolution to 
detect blends of Fe inner-shell absorption transitions in the 
ionized outflows from AGN for the first time. Unresolved ar- 
rays of inner-shell (2p-3d) transitions in moderately ionized 
PgO-15-i- jjj^yg been observed in the spectra of several type I 
A GN (see e .g. Behar et"al] ([2001') and references therein). 
iBehar et all (|200 calculate Fe UTAs for an ionization pa- 
rameter in the range log^ ~ 0.1 — 1.1. The UTA appears as 
a jagged, broad trough between ~ 16 — 17A. As ^ increases, 
the centroid of the UTA trough shifts towards shorter wave- 
lengths, the dominant contributions come from more highly 
ionized states of Fe and the Fe UTA broadens. 

Of the ten AGN in our sample that exhibit signatures 
of absorption due to photoionized material, only two AGN 
(NGC 4593 and NGC 5548) did not exhibit statistically sig- 
nificant evidence for a low ionization (Fe UTA) absorber 
component. Note that in the case of NGC 3516, we only 
found a low ionization component at statistical significance 
in the high fiux state of this AGN. Figure [5] shows the best- 
fit model of the UTA in each of the seven AGN with a warm 
absorber component where log^ < 1.1 (see Table [3} . The 
wavelength scale in Figure [S] corresponds to the outfiowing 
warm absorber rest-frame in each case. 



4.3 Soft X-ray emission in AGN 

As can be seen from Figs [2j [3] and |4l narrow emission lines 
were apparent in several of the AGN in our sample. The 
spectra of NGC 3783, NGC 4051 and NGC 5548 displayed 
the most prominent narrow lines. Narrow emission lines 
were less prominent in the spectra of Mkn 279, NGC 4593, 
Mkn 509 and MCG-6-30-15. The spectra of NGC 3227, 
NGC 3516 and NGC 7314 showed broad emission features 
around ~ 0.55 — 0.6 keV. However the broad emission fea- 
tures in the soft X-ray spectra of these three heavily ab- 
sorbed AGN were present in the respective background spec- 
tra and are therefore not intrinsic to the AGN. We tested 
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Figure 2. MEG observed photon spectrum compared to the best-fitting simple photoionized absorber plus emitter model (solid line) 
in the energy range 0.5-3.0 keV for NGC 3783 (see Tables |3] and |5]l . The data were binned at 0.02A for clarity. There are considerable 
gaps in the XSTAR database of Fe transitions so multiple Fe absorption features are not accounted for by the model. As can be seen, 
the absence of these lines from the model fit serves to alter the continuum fit in several places. There are also known problems with the 
fit. The continuum fit is poor around 2.0-2.5 keV (~ 5-6 A) in general due to calibration uncertainties. The continuum excess relative to 
the model around ~ 0.8 keV (~ 13-15 A) was also found by N03, using a different photoionization code (ION) to model the data. 
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Figure 3. MEG observed photon spectra compared to the best— fitting simple photoionized absorber plus emitter models (solid line) for 
five of the AGN in our sample (see Tables [3l and [5t . The data were binned at O.OSA for clarity in the plots (but 0.02A for fitting). The 
plot for NGC 3516 is that for the higher flux state, which is the one that showed evidence for a Fe UTA. Note that the emission feature 
at ~ 0.55 keV (observed frame) in NGC 3516 is duo to the background, not the source (see text). 
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Table 2. Best-fitting simple broken power-law continuum parameters (including neutral absorption with column 

density Nj{, covering a fraction, CF, of the source), and observed 0.5-2 keV fluxes. The broken power-law has 
photon indices Fi and F2, and a break energy Eb- 

Source cold Nh CF Fj Eg F2 Flux (0.5-2.0keV) 

(102"cm-2) (keV) (10^12 erg cm-^s-i) 

F9 7l 2.10t«;°? 1.06t°;°^ 1.78lg-°2 ^ 

3C 120 <1 l-79lo;lJ 0.99to^g 1.68 ±0.03 16.8 

NGC 3227 188 ±24 0.84l°;!^^ ' ' l-32l!5:i8 ^-^ 

NGC 3516 (low) 247t^^ 0.82^0.03 ^-^^-oaI 2.0 

NGC 3516 (high) 213^5^ 0.79t° ji[^ ^-^"^t^.H 3.7 

NGC 3783 102^2 0.75 ± 0.01 3.07lo;og 1.23 ± 0.01 l-92lo;g3 11.8 

NGC 4051 <1 3.08lo;j7 l-23to °g 1.77 ±0.05 12.6 

Mkn 766 2±1 3.10l°:28 0-9lio.04 l-SS^aM 15.0 

NGC 4593 3t_l 2.721^;^^ 0.96^0.02 1.80 ±0.04 20.9 

MCG-6-30-15 6± 2 3.821^;^^ ^■02^0°^'^ ^-^^-oTs 9-7 

IC 4329A 3213 3.661q;22 0.98 ± 0.03 l-83j;o;°3 77.9 

Mkn 279 <1 l-96lo:i6 0-89-0.09 1-56 ±0.03 5.6 

NGC 5548 <1 2.80lo:l3 0.95 ± 0.03 l-46lo;°3 ^0.0 

Mkn 509 <1 2.501q:o9 0-96to.06 1-63 ±0.03 23.1 

NGC 7314 107 ±5 0.97 ± 0.01 2.00to;og ^-^ 

Akn 564 4^1 3.16j;o;i5 l-43to o5 2.48 ± 0.05 44.9 




NGC 3783 



NGC 4051 



MCG-6-30-15 




Absorber Rest Wavelength (A) 



Figure 5. Plots of the model Fe UTAs produced by warm ab- 
sorber components from Table |3] with best-fitting ionization pa- 
rameters in the range log^ ~ 0.0 — 1.0 (sec text). 



for O VII and Ne ix triplet line emission and O vill Lya 
and Ne x Lya narrow lines in the spectra of the ten AGN 
listed above. We followed the fit procedure in ii4.1l except 
that the Gaussian model component was not inverted. In 
a given AGN, we established the weighted-mean offset ve- 
locity of the centroids of emission lines detected at > 90% 



significance in the spectrum. We then fitted XSTAR emis- 
sion line model grids at the mean offset velocity to those 
AGN spectra displaying multiple (> 2) emission signatures 
at > 99% significance. 

From the definition of ^, the emittin g material lies at a 
distance i? ~ 3 x 10~^-y/ Lion, 44/ (^100 ^^e) pc from the radia- 
tion source, where Lion, 44 is the 1-1000 Ryd ionizing radia- 
tion from the AGN (in units of lO** erg s~^), ^loo = C/100 
ergs cm~^s~^) and Ue is the electron density. For a volume 
filling factor of unity, the emitting material has a thickness 
AR ~ 1.2Nh /ne where Ue ~ 1.2n_ff and uh is the density of 
Hydrogen nuclei (assuming He/H ~ 0.1 and He contributes 
two electrons). Therefore, if we assume that AR < R, then 
rie > 1.7iV|iCioo/L44cm"^ where iV2i = Nh /W^^cm-^. The 
best-fit values of Nh and ^ from the XSTAR emission model 
fits then allow us to establish a lower limit on the electron 
density in the emitting material. It is important to note how- 
ever, that the emitting material is not necessarily the same 
as the absorbing material and therefore may not have the 
same physical parameters or location as the outflow respon- 
sible for absorption. 

The three AGN spectra with the most prominent 
narrow emission lines were NGC 4051, NGC 5548 and 
NGC 3783. In the spectrum of NGC 4051, the [O vil] and 
[Ne ix] forbidden lines were the most prominent, with addi- 
tional emission due to O vill Lya and Ne x Lya and weak 
O VII (i) emission. There was no apparent O vil (r) emis- 
sion, indicating that the emitting material in this source is 
photoionized rather than coUisionally ionized. Of course it 
is also possible that O vil (r) emission may have been re- 
absorbed by O VII (r) line absorption, which depends on 
the relative velocity shifts of the emitting and absorbing 
gases as well as the blending effects of the instrument res- 
olution and the spectral binning. The forbidden lines were 
marginally blueshifted from systemic velocity and the other 
(less prominent) emission lines were consistent with emis- 
sion at systemic velocity. The weighted mean velocity offset 
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Table 3. Detailed spectral fitting results for Photoionization models . Column 2 shows the best-fit cold 
absorption. Column 5 shows the mean offset velocity (negative=blueshift w.r.t. systemic) for the warm 
absorber components. Column 6 shows rate of mass outflow in terms of C, the product of the filling factor 
& the covering factor. Column 7 shows the fit-statistic (which includes emission model components from 
Table [Sj. Column 8 shows the increase in fit-statistic when the absorber component (3 dof) is removed 
from the model fit and the model is re-fitted . The true 90% confidence upper or lower limit could not 
be determined in these cases. The effective range of the offset velocity should be taken to lie between 
the best- fit value and the given upper or lower limit. 



Source 


Nh 


^ wabs 


log? 


'^wabs 


Mout/C 


C-stat 


A C 




(lO^^cm-^) 


(lO^icm-^) 


ergs cm s~^ 


km 


Mg/yr 






NGC 3516 (low) 


250 ± 20 


3±1 


2 4+"'2 
—0.1 


-910li« 


1.1 


1363 


60 


high flux 


2101^0 




2.5 ±0.1 


-5851^^ 


2.2 


1414 


131 






2±1 




-1685ll?5 


210 




18 


NGC 3783 


1021^ 


20+^ 


2.9 ±0.1 


-505 ± 15 


2.3 


9091 


689 






6±1 


2.1 ±0.1 


-515 ± 15 


14 




5070 






n + 0.2 
^-0.1 


0.4 ± 0.1 


_ . r+30 

-545+^^ 


720 




1039 






c+0.4 
^-1 




-il45_3(, 


0. 






NGC 4051 


< 1 


n 1 +0 1 

o-i-o:i 

9+?,1 
—0.3 


1.0 ± 0.3 

3.8 lb 0.1 


-520tgo 

-600+^? 
-2230tgo 


2.5 
0.1 
< 0.1 


1296 


29 
76 
41 




2 it 1 


— 0.5 


< 0.6 
3.1+?,'? 

— 0.1 


7t; + 25 
'O_70 

2^—145 


310 
< 0.1 


1289 


g 
6 








2.0 ±0.1 


~25_25 


0.4 




41 


NGC 4593 




2±1 

4+8 


2 4+°-^ 
^■^-0.2 
■1 O+0.2 
'^■■'-0.4 


yo_io 
-95^30 


0.3 
< 0.1 


1393 


25 
43 


MCG-6-30-15 


6 ± 2 


4±1 


0.2 ±0.1 




9.1 


1673 


226 






3± 1 


2.1±0.1 


-1 555+*° 


< 0.1 




240 








n 7+0.1 
■^■'-0.3 


0.2 




62 


IC 4329a 


32lt 


9 q+0.2 
^"^-0.3 


0.2 ±0.1 


-lOOlgs 


750 


1448 


162 






1.4 ±0.3 


2.2 ±0.1 


-100^20 


7.8 




175 


NGC 5548 


< 1 


0.6 ±0.2 


2.2 ±0.2 




89 


1065 


126 






50 ±40 


0.»_g 2 


_830+270 


3.3 




105 


Mkn 509 


< 1 


0.5 ±0.2 


9 O+0.2 
^■■'-0.1 


-140l:o 


14 


1137 


36 






0.1 ±0.1 




-1401^00 


1100 




27 


Akn 564 




0.1 ±0.1 


< 0.4 


-1401^0 


930 


1354 


18 






2+'^-^ 
^■^-0.1 


2.6 ±0.2 


-140+^5 


5.4 




131 



of the emission lines in this source (— 160 ± 75 km s ^) was 
apparently blueshifted. 

NGC 5548 also exhibited strong [O vil] and [Ne ix] for- 
bidden emission lines, with relatively strong O vill Lya and 
O VII (i) emission. O vil (r) line emission was weak in this 
source, indicating that the emitter is predominantly pho- 
toionized rather than coUisionally ionized. The He-like Ne 
and O triplet emission lines in NGC 5448 were marginally 
blueshifted ([O vil] lay at an offset velocity of —265^215 km 
s~^), but the O VIII Lya emission was marginally redshifted. 
The weighted mean offset velocity of the emitting material 
(—115 ± 125 km s~^) was consistent with emission at sys- 
temic velocity. NGC 3783 exhibited very strong O vill Lya 
, [Ne ix] and O vil (r) emission, the latter indicating that 
collis ional ionization may be important in the emitting ma- 
terial lPorguet et al.l (|200ll ). although photoexcitation due to 
high levels of UV flux could also account for O vil (r) emis- 



sion |Kinkhabwal^^e£^i] l|2002l ). Also present were emission 
features due to [O vii] and O vii (i) . The He-like forbid- 
den emission lines in NGC 3783 were marginally blueshifted, 
but the other emission lines were marginally redshifted. The 
weighted mean offset velocity of the emitting material was 
consistent with emission at systemic velocity. 

The spectra of Mkn 279, NGC4593, Mkn 509 and 
Mkn 766 each exhibited several emission features, most 
prominently O vill Lya and/or [O vil] , but they did not 
exhibit > 2 emission features at > 99% significance. Of the 
ten AGN, we therefore only fit the spectra of NGC 3783, 
NGC 4051 and NGC 5548 with XSTAR models of emis- 
sion. The best-fit results are given in Table [5] The re- 
quirement that AR/R < 1 yields lower limits of > 
3.7 X 10'', 2.3 X 10^ 2.8 x 10^ cm"^ for NGC 4051, NGC 5548 
and NGC 3783 respectively. 
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Table 4. Results from the literature for the spectra of the type I AGN in our study. For brevity, we only include the larger 
studies that parameterized Afj:/,log^ and/or the velocity of the warm absorber. We abbreviated X-ray instruments as: 
'KI{.=XMM-Newton RGS, XE=XMM-Newton EPIC-pn, CL=Chandra LETGS , CH=Chandra HETGS , By.=BeppoSAX , 
Ii\=RXTE . Gaps indicate value of log^ not specified in the literature. In these cases, either no photoionization modelling 
was carried out (Mkn 766, NGC 4051, NGC 3516) or a diff'erent form of the ionization parameter was specified (U or 
Uox)^ with insufficient information to translate to values of logg. ° Includes a model continuum component representing 
relativistically broadened O vm Lyo or a radiative recombination continuum from material close to a Kerr BH. Continuum 
model included a column density of 3 X lO^^cm"^ of Fe02. Continuum model included a blackbody with temperature 
of kT = 473 eV and a neutral absorber with a column density of ~ 3 X lO^^cm^^. Used z = 0.01676 from optical 
measurements rather than z = 0.01717 from 21cm H I measurements. . 



Source 


Instrument 


-^wabs 

{102lcm-2) 


log? 

(erg cm s"-*-) 




velocity 
km s~^ 


Reference 


NGC 3516 


XR 
CL 

XE,CL,BX 


~ 8 
~ 8 
16 
3-10 






1100 

1100 


Turner et al. ("2003") 
Netzer et al. (2002) 
Turner et al. (2005") 


NGC 3783 


CH 


20.0+^14^ 

10 0+^-^ 

7 9+2-1 
'■^-1.9 




[-400,-600], [ 
[-400,-600], [ 
[-400,-600], [ 


-1000,-1300] 
-1000,-1300] 
-1000,-1300] 


Netzer et al. (20031 


NGC 4051 


CH 
XE 

XE 


~ 1 
~ 0.1 
- 200 
^ 6 
~ 2 
[0.03,1] 


- 3.8 
2.7 ±0.1 
1.4 ±0.1 

[0,2.8] 




-2340 ± 130 
-600 ± 130 
~ -6500 
± 200 
0±200 
400 ± 100 


CoUinee et al. (2001) 
Pounds et al. (2004) 

Oele et al. (2004) 


Mkn 766 


CH 








0± 160 


Sako et al. (2003)'" 


NGC 4593 


CH,RX 
XR,CL 


5 4+1-5 

1.6 ±0.4 
0.1 ±0.05 


~ 2.5 
2.6 ±0.1 
0.5 ±0.3 




-135 ±40 
-400 ± 120 
-400 ± 120 


McKcrnan et al. (2003a) 
_StccnbruKge^^t_aL (2003b) 


MCG-6-30-15 


CH 
XR 


~ 32 
~ 5 
~ 2 
~ 2 


~ 2.5 
~ 0.7 
2-3 
0.5 - 2 




~ 
~ 

-1900 ± 140 
-150 ± 130 


Lee et al. (2001)'' 
Sako et al. (2003)'" 


IC 4329a 


XR 


~ 1.3 
~ 0.3 
6.6 ±0.4 
2.0 ±0.5 


-1.4 ±0.1 
0.6 ±0.1 
~ 1.9 
2.7 ±0.1 




~ 
[-100,-300] 
[-200,0] 
[-140,+180] 


Steenbruege et al. (2005a)'^ 


NGC 5548 


CL,CH 


2.5 ±0.5 


2.3 ±0.1 




-530 


Steenbruese et al. (2005b)'* 






1.0 ±0.4 
6+°-^ 

U.D_o 4 


1.9 ±0.1 
-0.2 ±0.2 




-530 
-530 




Mkn 509 


CH 


^•1-0.5 






-200 ± 100 


Yaaoob et al. (2003) 


Akn 564 


CH 


~ 1 

~ 1 


~ 2 
~ 1 




200 

~ -200 


Matsumoto et al. (2004) 






Table 5. Properties of Emitters . " We used Tie = lO^cm"'^ as 
a default electron density since the data could not discriminate 
between emission model grids in the range lO^^Hcm"''. 






Source 


JVem log 5em 

(10^"'^cm~2) (ergs cm/s) 


(cm-3) 


(km /s) 






NGC 3783 
NGC 4051 
NGC 5548 


13- '-13. 6 


1 61+°-^'^ 

1-°1-0.16 
1 40+°-^^ 


> 280 
> 3.7 X 10^* 

> 230 


-5 ±45 
-160 ±75 
-115 ± 125 



12 B. McKernan, T. Yaqoob & C. S. Reynolds 




10 20 

Obsen.^ed Wavelength (A) 



Figure 4. As for Fig. 3 for four more of the AGN in our sample. 



5 PROPERTIES OF THE WARM ABSORBERS 
IN TYPE I AGN 

The results of the photoionization modeUing of the AGN 
spectra are summarized in Table |3l Figures [2j [3] and |4] 
show the best-fit photoionization models (from Table |3} su- 



perimposed on the AGN spectra. Our uniform analysis of the 
AGN sample has yielded warm absorber parameters {Nh, 
log ^, velocity) that agree reasonably well with those in the 
literature (often from much more detailed analyses). 

Our results provide the first overview of high resolution 
soft X-ray grating data from Chandra observations of type I 
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AGN. In summary, we found that nine of the fifteen AGN 
in the sample have intrinsic continua that are more complex 
than a simple power law in the 0.5-5.0 keV band (modelled 
here as a broken power law). Of the remaining six AGN 
spectra, the continua of three (F9, 3C 120 and Mkn 279) are 
well described by simple power law models and the continua 
of the three remaining AGN (NGC 3227, NGC 3516 and 
NGC 7314) are well described by heavily absorbed power law 
models. Note that previous, higher S/N data of the heavily 
absorbed AGN (NGC 3227 & NGC 7314), evidence for pho- 
toio nized absorption has in some cases b een detected (see 
e.g. lPtak et al] l| 19941 ) ICeorge et all (|l998D ). The nine AGN 
spectra with complex intrinsic continua exhibit signatures 
of an ionized absorber as does the spectrum of NGC 3516, 
in spite of being heavily absorbed. The absorbers appear 
to be photoionized and outflowing, with velocities in the 
range ~ — 2000 km s^^ (similar to the speed of the So- 
lar wind). The column density of the warm absorbing gas 
is ~ 10^''~^^cm"^. Nine of the ten AGN exhibiting warm 
absorption are best-fit by multiple ionization components 
and three of the ten AGN require multiple kinematic com- 
ponents. The warm absorbing gas in our AGN sample has 
a wide range of ionization states, spanning roughly four or- 
ders of magnitude (^ ~ lO'^ — 10*) and up to three orders 
of magnitude in the same source. Of the ten AGN spec- 
tra that exhibit warm absorption, a simple fit to the 0.5-0.7 
keV and 2.5-5.0 keV energy regions reveals that the intrinsic 
continuum of two AGN (NGC 4051 and NCG 5548) exhibits 
a statistically significant upturn in the intrinsic soft X-ray 
continuum relative to a hard X-ray power law. Our simple 
test shows that MCG —6—30—15 does not require a soft ex- 
cess. The remaining AGN that exhibit warm absorption are 
either heavily absorbed in soft X-rays or exhibit spectral 
complexity that cannot be accounted for by such a naive 
test of the data. Nine of the ten AGN that exhibit warm 
absorption also exhibit Fe UTAs, indicative of an absorber 
component that could (depending on absorber density and 
geometry) carry most of the mass in outflow, as we show 
below. 

Fig. |6] compares the warm absorber component param- 
eters from the different AGN (see Tables [3] & [S} . Fig. EI a) 
shows the relationship between warm absorber outflow ve- 
locity and the corresponding warm absorber column density. 
Clearly the warm absorber components in all ten AGN are 
outflowing and with velocities that span around two orders 
of magnitude. There appears to be a gap in the outflow ve- 
locities in our sample between ~ 300 — 500 km s~^, the 
origin of which is not clear. The outflow components with 
velocities below this gap tend to be associated with lower 
column densities than those with with velocities above the 
gap. Fig. [nib) shows that the ionization parameter of the 
warm absorber is independent of velocity and can vary by 
up to three orders of magnitude in the same source, sug- 
gesting that the outflowing absorber components consist of 
gas lying at different poin ts on the heating/cooling curve 
(e.g. Krolik fc Krisj l|2001^ '). Furthermore Fig. [gib) also in- 
dicates that there is a separate population of low ionization 
state absorber components at ^ < 10. These low ionization 
parameter absorber components are responsible for the Fe 
UTAs in the AGN spectra. One caveat is that there is no 
bias in principle in the X-ray band for absorbers with ion- 



(a) 



■XX- 



V (km s ) 



Figure 6. Velocity of Warm Absorber Components (a) vs Col- 
umn Density {Njj) and (b) vs ^ of Warm Absorber Component 
for the AGN in our sample (from Table [Sj . In panel (a) , we plot 
the total column density of all the warm absorber components 
around a particular velocity. Thus, e.g. the three low velocity 
components from NGC 3783 in Table [3l become a single velocity 
component between -490 and -565 km s~^ in panel (a). In panel 
(b), we plot all the individual warm absorber components to il- 
lustrate the fact that the ionization parameter can vary by up to 
three orders of magnitude in the same source. In panel (b), the 
low ionization components (^ < 10.) correspond to the Fe UTAs. 



ization states 10 < ^ < 100, however Chandra HETGS is 
less sensitive than other instruments to this band. 

The warm absorber outflow apparent in ten of the fif- 
teen AGN in our sample must clearly carry mass away from 
the central SBH. The rate of mass loss from an AGN outflow 
is given by 



Mout 



m„nev4:TTr 



(1) 



where AQ is the solid angle subtended by the absorber 
at the ionizing source, C is the volume filling factor, y is 
the mean atomic mass per Hydrogen atom {y ~ 1.3 for so- 
lar abundances), nip is the proton mass, is the electron 
density of the absorber, r is the distance of the absorber 
from the radiation source, and v is the outflow velocity. We 
believe that the filling factor, C, cannot be reliably con- 
strained by current observations of AGN (see e.g. discussion 
in H5.ip . Blustin et al. (2005) estimated C for a sample of 
sources based on various assumptions, but those same as- 
sumptions led Blustin et al. (2005) to derive a maximum 
distances between the absorber and radiation source that 
was less than the minimum distance for no less than five 
sources. On the other hand, it is possible to argue that the 
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covering factor (Af2/47r) is approximately equal to the frac- 
tion of type I AGN in which a warm absorber is detected, 
or ~ 0.5 — 0.7. However, this argument breaks down if AGN 
winds are accelerating and/or "bend", rather than constant 
velocity, non-accelerating outflows. Even so, this still leaves 
C (and therefore the mass outflow rate) unconstrained. From 
the definition of the ionization parameter, = Lion/ne(,, 
therefore 



V 47r / f 100 



(2) 



where Lion,AA = Lion/ 10*"*, Cloo = C/100 and nsoo = 
f/SOO where v is the outflow velocity, and x is the number 
of Hydrogen atoms per free electron. With x — 9/11 and 
y — 1.3, the product xy is approximately unity, or ~ 1.06. 
Equation [2] suggests that all else being equal, warm absorb- 
ing gas with a low value of ionization parameter should ac- 
count for the largest mass outflows from AGN. From Ta- 
ble [3] we flnd mass outflow rates in the range ~ 10"^ -10^ X 
(An/47r)CM0yr"\ Therefore, even for small filling factors, 
the outflow rate can be comparable to the expected accretion 
rate onto the central supermassive black hole. The accretion 
rate onto the black hole can is Maccretion = L^ao\M(? , where 
Lboi is the bolometric luminosity of the AGN and rj is the 
accretion efficiency of the black hole. If Lboi = 10XioLion,44, 
where Xio is a parameter of order unity, and r] = O.I770.1, 
then A^accrction ~ Li on,44-'fio/(5.7 X rio.i)MQ/yr and 



~ 94 ( ^ ) ( ^ ) ( ^ ) r^o.iC 

Maccretion V 47r / V AlQ / V ^lOO 



(3) 



The value of C and the covering factor of the absorber 
is critical to understanding the processes underpinning the 
warm absorption phenomenon. If W500, ^10, ??o.i ~ 1 and if 
(An/47r) ~ 0.5 typically, then the filling factor must be very 
small (< 0.02) for the mass outflow rate to be comparable 
to or less than the accre tion rate (a point also made by 
ISteenbrugge et all (|2005bl )). 

The warm absorber components found in Seyfert 1 
galaxies have column densities that are too small to account 
for the EWs of the narrow, core, Fe K line emission. How- 
ever, for our sample, we investigated whether the core Fe K 
line emission EW correlate with the column densities of the 
warm absorber components in case of a secondary effect. 
For example, if the warm absorber material is an outflow 
of material from the putative obscuring torus, and if the 
core Fe K line emission originates in the torus itself, one 
might expect a correlation between the Fe K line EW and 
the warm absorber column density if thicker torii produce 
thicker winds. We also investigated whether the Fe Ka line 
in Type I AGN 'knows' about the warm absorber by com- 
pared the total column density of the warm absorber in an 
AGN with the EW of the corresponding core Fe Ka line 
l|Yaaoob fc Padmanabh an 2004) (Figure [T]). 

For our sample. Figure [7] shows the EW measurements 
(from Yaqoob & Padmanabhan 2004) versus warm absorber 
column density. A weak correlation is permitted based on 
Figure [T] so the Fe Ka line might 'know' about a large 
warm absorber column density. In this case we calculate 
the Pearson correlation coefficient to be 0.88 (for N=9), 
which is significant at the ~ 99% confidence level. How- 
ever, we caution that the EWs of individual Fe Ka lines 
can vary by a factor of 2-3 as the continuum varies whilst 



50 100 
EW FeK (eV) 



Figure 7. EW of the Fe Ka line in each AGN versus the total 
warm absorber column density (from Table [5|l. Fe Ko: measure- 
ments were taken from I Yaqoob fc PadmanabhanI ll2004h . A weak 
correlation between total column density of absorber and the EW 
of the Fe Kq line is permitted, however Fe Kq line EWs can vary 
significantly in individual AGN, so a correlation is unlikely. 



the line intensity does not, and this would destroy any hint 
of correlation in Figure [T] which may therefore be spuri- 
ous. In particular, we note that no warm absorber is de- 
tected in F9, yet it has a promin ent, strong Fe Ka line core 



a promm e 

dYaaoob fc PadmanabhanI 120041 '). 



5.1 Comparison with other AGN sample studies 

iBlustin et all l|2005l ) collated from the literature the results 
of some of the X-ray spectral observations of 23 AGN using 
Chandra and XMM-Newton and attempt to infer the general 
properties of Seyfert warm absorbers from these results. The 
methodology of the iBlustin et all (|2005l ) study is quite dif- 
ferent from the present work. However, the strengths of our 
study (uniform analysis and a single instrument) may also 
be weaknesses. Our selection criteria for warm absorption 
signatures may be too conservative, our uniform modelling 
of the continuum may ignore individual source peculiarities 
and our band-pass may miss e.g. very low ionization ab- 
sorber components. A study such as the one carried out by 



Blustin et al 



Blustin et al 



2005 ) is a useful comparison for these reasons. 



20051 ) conclude that most warm absorbers are 



most likely to originate in outfiows from the dusty torus, 
since they estimate that the minimum distance of the warm 
absorbers from the central radia tion source mostly c luster 
around the distance of the torus. IBlustin et all (|2005l l con- 
clude that the kinetic luminosity of these outfiows is < 1% 
of the AGN bolometric luminosity and the observed soft X- 
ray absorbing ionization phases fill < 10% of the available 

volume. 

Some of the conclusions reached bv lBlustin et al.1 (|2005l ) 
are quite different from ours. We make no claim for the flU- 
ing factors of the outflows or indeed, their origin. First, we 
believe there are currently insufficient data and constraints 
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Table 6. Black hole masses for AGN with warm absorbers . Column 2 lists 
the 0.5-2.0 keV luminosi ty for each sourc e, corr ected for absorption. Refere nces: 
(1) [P eterson ct al.' ('2004'). (2) iBotte et al.] ||2005|) . (3) iMcHardv et al.l |200l). (4) 
iBian fc Zhao (2003). 



Source 


0.5-2.0 keV 
luminosity 

(lO^^org/s) 


Lion 

(lO^^crg/s) 


Mbh 
(10«Mq) 


Lion/ L^dd 


Reference 


NGC 3516 (low) 


0.16 


0.15 


43 


0.003 


1 


NGC 3516 (high) 


0.38 


0.36 




0.007 


1 


NGC 3783 


1.25 


1.00 


30 


0.003 


1 


NGC 4051 


0.02 


0.02 


2.0 


0.006 


1 


Mkn 766 


0.66 


0.47 


0.6 


0.572 


2 


NGC 4593 


0.50 


0.29 


5.4 


0.042 


1 


MCG-6-30-15 


0.32 


0.16 


4.5 


0.027 


3 


IC 4329A 


0.76 


3.40 


7.5 


0.349 


1 


NGC 5548 


6.98 


8.33 


96 


0.067 


1 


Mkn 509 


6.55 


5.23 


100 


0.040 


1 


Akn 564 


6.54 


4.69 


1.2 


3.000 


4 



to reliably estimate the filling factor, preventing the deriva- 
tion of absolute ma ss outflow rates. Second, we note that 
iBlustin et all (|2005l ) calculate a minimum distance between 
the central radiation source and the warm absorber based 
only on the assumption that the outflow velocity exce e ds the 
escape velocity. However, very recently I Miller et al.l (|2006l ) 
observed an X-ray absorbing outflow from the stellar-mass 
black hole binary GRO J1655-40 with a radial velocity far- 
less than the e scape velocity at its location. Moreover, the 
assumptions of lBlustin et al.l ()2005l ) lead to a maximum dis- 
tance of the warm absorber from the radiation source that is 
less than the minimum distance for no less than five AGN. 
Robust information on the distance between the warm ab- 
sorber and the central source can only come from var iabil- 
ity studies of the absorber. Third, iBlustin et al.l (|2005[ ) con- 
clude that Seyfert warm absorbers are probably not telling 
us anything fundamental about the energetics or structure of 
the central engine. However, if for example, warm absorbers 
originate in a disk wind, the properties of the disk wind 
are likely to correlate with those of the disk (magnetic flux 
strength, disk temperature gradient, velocity of disk etc). Fi- 
nally, some parts of t he warm absorb ing outflows from AGN 
discussed here and bv IBlustin et al. I ((2005. ) may not actually 
he observable because they are fully ionized. Such fully ion- 
ized outflows may be the fastest component of the outflow 
since it may be closest to the source of ionizing radiation. 



6 DISCUSSION 

At present there are two distinct th eories of the origin of 
the warm absorber. On the one hand, lElvid (|2000f ) proposes 
a unified scheme for quasars that includes a wind that rises 
vertically from a narrow range of radii on the inner accre- 
tion disk and is then bent outwards by a radial radiation 
force to produce a funnel-shaped thin shell outfiow. This 
model specifically excludes the dusty molecular torus as a 
characteristic of AGN and pro poses that t he Fe Ka emission 
line is produced in the funnel. lElvid l|2000l ) predicts that the 
range of line widths will be similar to the broad absorption 
line 'detatchment velocities' (~ — 5000 km s^^ vertically, 
so one might expect a narrower line width than this when 



looking down the funnel). Of the forty-eight absorption fea- 
tures used to constrain the warm absorber components in 
this study, thirty-seven have FWHM < 2000 km s"^ and 
only four allow for FHWM > 3000 km s"^ Statistically, a 
weak correlation is allowed between the EW of the FeKa line 
and the total column density of the warm absorber. How- 
ever, we note that in our sample, the source with the least 
evidence for any warm absorption (F9) exhibits the largest 
FeKa EW, which is not consistent with expectations from 
the Elvis model. 

On the other hand, iKroUk fc Kr^ (|200ll ) propose that 
the warm absorber originates in a photoionized evaporation 
from the inner edge of the putative obscuring torus believed 
to surround the AGN central engine. In this model, the warm 
absorber is a multi-temperature wind, with the different out- 
flow components having values of H (pressure form of ion- 
ization parameter) that cover a relatively narrow range, as 
a result of the different phases of gas co-existing in pressure 
equilibrium. From our photoionization model fits, we find 
that the temperatures of the warm absorber components lie 
in the range T ~ 10*~^ K for the sample, yet logH lies in 
the range ~ — 1, with the highest ionization components 
generally having the largest values of T and S. Thus, our 
general constr aints seem to agree qu ite well with the model 
predictions of iKrolik fc Kriss (|200ll ). An additional attrac- 
tion of this model is that the interaction of the disk wind 
with the dusty molecular torus in this model can also nat- 
urally explain dusty warm absorption, which may be the 
cause of additional spectral complexity of some of the AGN 
in this sample (see 

One important point to note is that some of the warm 
absorber components discussed above may actually be due 
in part to hot local gas at cz~0. For example the low ve- 
locity warm absorber component in NGC 4051 at (~ —600 
km s~^) is kinematically very close to an absorption signa- 
ture at cz~0 since NGC 4051 is cosmologically redshifted by 
cz = 726 ± 15 km s~^ from z~0. Other AGN in this sample 
yield tantalizing hints of abso rption due to hot local gas (see 
iMcKernan et all (|2004l . [2OO5I ) and references therein for fur- 
ther details). Of course, the absorption signatures of some 
of the warm absorbing outfiows could be mimicked by hot 
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gas at intermediate redshift. However, the expected column 
density of most filame nts of warm/hot intergalactic medium 
l|Cen fc Ostrikerlll999l l is far lower than that observed by us. 

The soft X-ray spectra of type 1 AGN observed with 
the high resolution Chandra and XMM-Newton observa- 
tories prompt many intruiging questions. We believe that 
these questions can be answered only by the careful uni- 
form analysis of a sample of AGN. Our uniform analysis 
has its limitations since we use spectra from only one in- 
strument with a limited bandpass, and our methodology is 
sufficiently conservative that we may miss less signficant ab- 
sorber signatures. Nevertheless, we have carried out a first 
uniform analysis on a small sample of AGN observed with 
Chandra HETGS and we have established reasonably well 
constrained parameter ranges for the warm absorbing out- 
flows in type I AGN. Interestingly, we found that mass loss 
resulting from the warm absorber outflow can be high, com- 
parable to or greater than the expected accretion rate onto 
the central supermassive black hole. Most of the outflow- 
ing mass could be carried by low ionization state outflows 
(depending on absorber density and geometry), which are 
best studied in the X-ray band via their Fe UTA spectral 
imprint. Low rates of mass outflow from AGN and higher 
velocity outflows may be associated with lower values of the 
AGN Eddington ratio, but the rate of mass outflow does not 
correlate with the mass of the central black hole. A weak 
contribution to the core of the narrow Fe Kq line from high 
column density warm absorbers seems unlikely. There is gap 
in the outflow velocities of the warm absorber components 
in the range ~ 300 — 500 km s~^ that is puzzling. All these 
results so far provide tantalizing hints at global patterns of 
type I AGN behaviour and merit further investigation using 
a larger sample of high spectral resolution AGN spectra. 



ACKNOWLEDGMENTS 

BM gratefully acknowledge support from NSF grant 
AST0205990. TY acknowledges support from NASA 
through grant AR4-5009X issued by Chandra X-ray Ob- 
servatory Center, operated by the SAO for and on be- 
half of NASA under contract NAS8-39073. We made use 
of the HEASARC on-line data archive services, supported 
by NASA/GSFC and also of the NASA/IPAC Extragalactic 
Database (NED), operated by the Jet Propulsion Labora- 
tory, CalTech, under contract with NASA. Thanks to the 
Chandra instrument and operations teams for making the 
observations possible. Thanks to Tim Kallman for numer- 
ous useful discussions on XSTAR. Thanks to the anonymous 
referee for very detailed and useful comments that helped 
improve and shorten this paper. Thanks to Fabrizio Nicastro 
for discussing his results from an unpublished XMM-Newton 
RGS observation of NGC 4051. 



REFERENCES 

Arav, N., Korista, K. T. & de Kool, M. 2002, ApJ, 566, 
699 

Behar, E., Sako, M., & Kahn, S. M. 2001, ApJ, 563, 497 
Behar, E. & Netzer, H., 2002, ApJ, 570, 165 
Bian, W. & Zhao, Y. 2003, MNRAS, 343, 164 



Blustin, A. J., Branduari-Raymont, G., Behar, E., Kaastra, 

J. S., Kahn, S. M., Page, M. J., Sako, M., Steenbrugge, 

K. C, 2002, A&AS, 392, 453 
Blustin, A. J., Page, M. J., Fuerst, S. V., Branduardi- 

Raymont, G. & Ashton, C. E., 2005, A&AS, 431, 111 
Botte, v., Ciroi, S., diMille, F., Rafanelli, P. & Romano, 

A., 2005, MNRAS, 356, 789 
Branduardi-Raymont, G., Sako, M., Kahn, S. M., 

Brinkman, A. C, Kaastra, J. S., & Page, M. J. 2001, 

A&AS, 365, 140 
Cash, W., 1976, A&AS, 52, 307 
Cen, R & Ostriker, J. P., 1999, ApJ, 514, 1 
Collinge, M. J., Brandt, W. N., Kaspi, S., Crenshaw, D. M., 

Elvis, M., Kraemer, S. B., Reynolds, C. S., Sambruna, R. 

M., Wills, B. J. 2001, ApJ, 557, 2 
Crenshaw, D. M., & Kraemer, S. B., 1999, ApJ, 521, 572 
Elvis, M., Wilkes, B. J., Lockman, F. J. 1989, AJ, 97, 777 
Elvis, M., 2000, ApJ, 545, 63 

Ferland, G. J., Martin, P. G., van Hoof, P. A., & Wein- 
garter, J. C. 2002, Workshop on X-ray Spectroscopy of 
AGN with Chandra and XMM-Newton, held at MPE 
Garching, December 3-6, 2001, Report 279, p. 103 

Gehrels, N. 1986, ApJ, 303, 336 

George, I. M., Turner, T. J., Netzer, H., Nandra, K., 

Mushotzky, R. F., Yaqoob, T. 1998, ApJS, 114, 73 
Halpern, J. P., 1984, ApJ, 281, 90 

Kaastra, J. S., Steenbrugge, K. C, Raassen, A. J. J., van 

der Meer R. L. J., Brinkman, A. C, Liedahl, D. A., Behar, 

E., & de Rosa, A. 2002, A&AS, 386, 427 
Kallman, T. 2003, private communication 
Kinkhabwala, A. et al. 2002, ApJ, 575, 732 
Krolik, J. H., & Kriss, G. A. 2001, ApJ, 561, 684 
Lee, J. C, Ogle, P. M., Canizares, C. R., Marshall, H. L., 

Schulz, N. S., Morales, R., Fabian, A. C, & Iwasawa, K. 

2001, ApJ, 554, L13 
Markert, T. H., Canizares, C. R., Dewey, D., McGuirk, M., 

Pak, C, & Shattenburg, M. L. 1995, Proc. SPIE, 2280, 168 
Mason, K. O. et al. 2003, ApJ, 582, 95 
Matsumoto, C, Leighley, K. M. & MarshaU, H. L. 2004, 

ApJ, 603, 456 

McHardy, L M., Gunn, K. F., Uttley, P., Goad, M. R., 

2005, MNRAS, 359, 1469 
McKernan, B., Yaqoob, T., George, LM., Turner, T. J. 

2003a, ApJ, 593, 142 
McKernan, B., Yaqoob, T., Mushotzky, R., George, LM., 

Turner, T. J. 2003b, ApJ, 598, L83 
McKernan, B., Yaqoob, T., Reynolds, C. S. 2004, ApJ, 617, 

232 

McKernan, B., Yaqoob, T., Reynolds, C. S. 2005, MNRAS, 
361, 1337 

Miller, J. M. et al. 2006, Nature (accepted), 

astro-ph/0605390 
Morales, R. M. & Fabian, A. C, 2002, MNRAS, 329, 209 
Murphy, E. M., Lockman, F. J., Laor, A., Elvis, M. 1996, 

ApJS, 105, 369 
Nandra, K. & Pounds, K. A. 1992, Nature, 359, 215 
Netzer, H., Chelouche, D., George, L M., Turner, T. J., 

Crenshaw, D. M., Kraemer, S. B. & Nandra K. 2002, ApJ, 

571, 256 

Netzer, H., et al. 2003, ApJ, 599, 933 
Ogle, P. M. et al. 2004, ApJ, 606, 151 
Peterson, B. M. et al. 2004, ApJ, 613, 682 



A Soft X-ray study of type I AGN observed with Chandra 



Porquet, D., Mewe, R., Dubau, J., Raasscri, A. J. J. & 

Kaastra, J. S. 2001, A&AS, 376, 1113 
Pounds, K. A., Reeves, J. N., Page, K. L., Edolson, R., 

Matt, G. & Perola, G. C. 2003, MNRAS, 341, 953 
Pounds, K. A., King, A. R., Reeves, J. N., Page, K. L. 2004, 

MNRAS, 350, 10 
Ptak, A. F., Yaqoob, T., Serlemitsos, P. J., Mushotzky, R., 

Otani, C. 1994, ApJ, 436, L31 
Reynolds, C. S. 1997, MNRAS, 286, 513 
Reynolds, C. S. & Nowak, M. A., 2003, Phys. Rep., 377, 

389 

Sako, M., et al. 2001a, A&AS, 365, L168 

Sako, M., ct al. 2003, Ap.J, 596, 114 

Steenbrugge, K. C, Kaastra, J. S., Blustin, A. J., 

Branduardi-Raymont, G., Sako, M., Bchar, E., Kahn, S. 

M., Paerels, F. B. S., Walter, R. 2003b, A&AS, 408, 921 
Steenbrugge, K. C, Kaastra, J. S., Sako, M., Branduardi- 
Raymont, G., Behar, E., Paerels, F. B. S., Blustin, A. J. 

and Kahn, S. 2005, A&A, 432, 453 
Steenbrugge, K. C. et al. 2005, A&A, 434, 569 
Turner, T. J., Nandra, K., George, I. M., Fabian, A. C, 

Pounds, K. A. 1993, ApJ, 419, 127 
Turner, T. J., Kracrncr, S. B., Mushotzky, R. F., George, 

I. M. & Gabel, J. R. 2003, ApJ, 594, 128 
Turner, T. J., Kraemer, S. B., George, I. M., Reeves, J. N., 

& BottorfT, M. C. 2005, ApJ, 618, 155 
Ward, M., Elvis, M., Fabbiano, N., Carleton, P., Willner, 

S. P., & Lawrence, A. 1987, ApJ, 315, 74 
Yaqoob, T., McKernan, B., Kraemer, S. B., Crenshaw, D. 

M., Gabel, J. R., George, I. M. & Turner, T. J., 2003, 

ApJ, 582, 105 
Yaqoob, T. & Padmanabhan, U., 2004, ApJ, 604, 63 



This paper has been typeset from a T^jX/ I^-T^iX file prepared 
by the author. 



